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Abstract 
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CMS detector, corresponding to an integrated luminosity of 19.7 fb^^. The h boson is 
assumed to be the standard model-like Higgs boson with a mass of 125 GeV. With no 
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section and the branching fraction of the A boson in the Zh channel. Results are also 
interpreted in the context of two Higgs doublet models. 
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1 Introduction 

The discovery of a scalar boson at the CERN LHC |Il]-|3l with properties in agreement with 
those predicted by the standard model (SM) raises the question of whether the Higgs sector 
consists of only one physical state, as expected in the SM, or whether additional bosons are 
also involved. 

An extension of the SM Higgs sector is provided in two Higgs doublet models (2HDM) |4J, 
which introduce a second scalar doublet in addition to the one from the SM. Different formu¬ 
lations of 2HDM predict different couplings of the two doublets to quarks and to leptons. In 
Type-I 2HDM, all fermions couple to only one Higgs doublet, while in Type-II, up- and down- 
type quarks couple to different doublets. One example of a Type-II 2HDM is the minimal 
supersymmetric standard model ||5l, despite that supersymmetry is not explicitly required in 
2HDM. The second Higgs doublet entails the presence of five physical states: two neutral, CP- 
even states, representing a light h and a heavy H boson; a neutral, CP-odd A boson; and two 
charged scalar H^ bosons. The lightest scalar h is assumed to be the boson observed at the LHC 
at a mass of 125 GeV IIHII. If the masses of the heavier bosons are at or below the TeV scale, 
they can be accessible at the LHC. Searches for this extended sector can be performed either by 
measuring the values of the couplings of the discovered h boson to other SM particles IITOl - IT^ . 
or via direct searches in final states disfavored by the SM fl^ . A way to probe this kind of 
new physics is therefore to search for bosons that decay into final states that contain an SM-like 
Higgs boson. 

This paper describes a search for a heavy pseudoscalar A boson that decays into a Z and an h 
boson, both on-shell, with the Z boson decaying into a pair of leptons {i being e or }i), and 
the h boson into bb. In most 2HDM formulations [4], the A boson is produced predominantly 
through gluon-gluon fusion and decays to on-shell Z and h bosons, provided that the mass of 
the A boson satisfies jma ^ ~ 216 GeV. This charmel is expected to be viable for wja 

smaller than twice the top quark mass (nit), where the decay A —)■ Zh is generally dominant, 
but, depending on model parameters, it can also be sensitive at larger values of nip^. The h —> bb 
decay has a large branching fraction for most of the parameter space in 2HDM IfTTl . A similar 
analysis has been recently published by ATLAS Ill4t . 

The analysis strategy is to reconstruct the Z, h, and A boson candidates from the visible decay 
products in the event. The signal would manifest itself as a peak in the four-body invariant 
mass spectrum over an expected SM continuum. Irreducible backgrounds correspond 

to Z boson production with two accompanying b quark jets, and tf events in the dileptonic fi¬ 
nal state. These backgrounds are evaluated and normalized directly using appropriate control 
regions in data. The h boson produced in association with a Z boson provides a contribution 
to the background, but it differs from signal because the mass does not contain a reso¬ 
nant peak. Signal sensitivity is improved by exploiting the known value of the h boson mass, 
using it to rescale the jet momenta to match the value expected for the dijet invariant mass. In 
addition, optimal signal efficiency and background rejection is achieved using a multivariate 
discriminator. Results are extracted through a two-dimensional (2D) fit to m^£bb srid the dis¬ 
criminator output; upper limits are presented on the product of the total cross section and the 
A ^ Zh, Z —> £i, and h ^ bb branching fractions for a pseudoscalar boson, and interpreted 
within the 2HDM. 
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3 Data and simulation 


2 CMS detector 

A detailed description of the CMS detector, together with a definition of the coordinate system 
and kinematic variables, can be found in Ref. IITSl . The central feature of the CMS apparatus 
is a superconducting solenoid of 6 m internal diameter, providing a magnetic field of 3.8 T. 
The field volume contains a silicon pixel and strip tracker, a homogeneous electromagnetic 
calorimeter (ECAL), and a sampling hadron calorimeter (HCAL), each composed of a barrel 
and two endcap sections. Muons are measured in gas-ionization detectors embedded in the 
steel flux-return yoke outside the solenoid. 

The silicon tracker measures charged particles within the pseudorapidity range of |//| < 2.5. 
For non-isolated particles with transverse momenta of 1 < px < lOGeVand \i]\ < 1.4, the track 
resolutions are typically 1.5% in px, and between 25-90 and 45-150 pm, respectively, in trans¬ 
verse and longitudinal impact parameters relative to the production vertex ||T6| . The ECAL 
consists of lead tungstate crystals that provide a coverage up to |p| < 3.0. The mass resolution 
for Z —>• e+e^ decays when both electrons are in the ECAL barrel is 1.6%, and is 2.6% when 
both electrons are in the endcaps. The HCAL has alternating layers of brass as absorber and 
plastic scintillators, and covers the range of |p| < 3.0, which is extended to |p| < 5.2 through 
forward calorimetry. Muons are measured in the range of | p | < 2.4, with detection planes 
made using three technologies: drift tubes, cathode-strip chambers, and resistive-plate cham¬ 
bers. Matching muons to tracks measured in the silicon tracker provides a px resolution of 
1.3-2.0% for muons with 20 < px < 100 GeV in the barrel and better than 6% in the endcaps. 


3 Data and simulation 

Data used for this analysis were collected using double-muon and double-electron triggers, 

with Px thresholds set to 17 and 8 GeV for the highest and next highest px lepton, respectively, 

and an isolation requirement used in the electron trigger to maintain an acceptable rate. The 

_1 

analyzed events correspond to an integrated luminosity of 19.7 ± 0.5 fb of pp collisions at 
v/s = 8TeV CZl. 

Signal samples and SM background processes Z+jets, W+jets, and tf+jets or a vector boson 
(tfV) are simulated using the MadGraph 5.1 IllSl Monte Garlo (MG) generator; multijets and 
dibosons {VV, with V, V being W or Z) are generated using LO PYTHIA 6.4 IfT^ and the 
GTEQ6L |[20l parton distribution functions (PDF). Single top quark production, and SM Higgs 
boson production in association with an electroweak vector boson (Vh) are generated using 
the next-to-leading-order (NLO) powheg 1.0 ll2Tl - (2^ MG generator and the MSTW2008NLO 
PDF Il24l . Parton showering and hadronization are performed with PYTHIA using the Z2* 
tune 1251 . The generated MG events, including additional pp interactions (pileup) occurring 
in the bunch crossing containing the high-px scatter, are processed through a full detector sim¬ 
ulation based on Geant4 Ii26l and reconstructed with the same algorithms as used for data. 

The pseudoscalar A boson is assumed to be produced via the gluon-gluon fusion process and to 
have a narrow width. The validity of the narrow-width approximation is discussed in Section]^ 
The branching fraction S (A —I Zh) is set to 100%, and similarly only the Z boson decays into 
electrons or muons, and the h boson decays into a pair of bb quarks, are considered. Other 
decay modes have negligible efficiency for passing the selections. The mass of the light Higgs 
boson is set to nij^ = 125 GeV, while the search for the A boson is performed in the mass range 
from 225 to 600 GeV, above which the efficiency to reconstruct two separate jets from h decay 
becomes too small because of its increased momentum. 
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4 Event reconstruction 

A global reconstruction of the event is achieved using a particle-flow (PF) technique, which 
reconstructs and identifies individual particles emerging from each collision using information 
from all CMS subdetectors ||2Zl|28l. 

Electron candidates are reconstructed for \ri\ < 2.5 by matching energy depositions in the 
ECAL with reconstructed tracks f29|. The identification relies on a multivariate technique that 
combines observables sensitive to the amount of bremsstrahlung along the electron trajectory, 
the compatibility of the measured position, direction and momentum of the reconstructed track 
in the irmer tracker, and the energy deposition reconstructed in the ECAL cluster. Additional 
requirements are imposed to remove electrons produced by photon conversions in the detector 
material. 

Muons are reconstructed within | < 2.4, combining information from both the silicon tracker 
and the outer muon spectrometer Il30l . requiring small energy depositions in the calorime¬ 
ters |3T]| . Muon candidates have to fulfill restrictive selection criteria based on the quality and 
the impact parameter of the track, as well as on the number of hits observed in the tracker and 
muon systems. 

An isolation variable is defined for each lepton through the scalar sum over the pj of all PE 
candidates, excluding the lepton, within a cone of AR = V {Arj)^ + < 0-4 around the 

lepton direction, where (p is the azimuthal angle in radians, subtracting the pileup contribu¬ 
tion 1^, and then dividing by the lepton pj. The lepton is rejected if the isolation exceeds 0.15 
for electrons and 0.12 for muons. 

Jets are formed from all particles, charged and neutral, reconstructed through the PE algorithm, 
and clustered with the anti-kj algorithm JSSj |3l| with a distance parameter of 0.5. Jet energy 
corrections are determined from dijet and Z/ 7 +jet events, and applied to both data and simula¬ 
tion [351 . The imbalance in transverse momentum is calculated as the negative of the vectorial 
sum of transverse momenta of all the PE candidates, and its magnitude is denoted as [361 . 

The combined secondary vertex (CSV) b tagging algorithm [371 is used to identify jets that 
originate from b quarks. This algorithm combines information from track impact parameters 
relative to the primary and secondary (displaced) vertices into a likelihood discriminant. Two 
standard working points are set on the discriminant, corresponding to restrictive (tight) and 
less-restrictive (loose) thresholds that provide on average 50% and 80% b tagging efficiencies, 
with respective misidentification probabilities of about 0 . 1 % and 10 % for light-flavor jets, and 
about 5% and 25% for c jets. The CSV distribution is corrected in simulation to take into account 
a difference at the percent level in algorithm performance for data and simulation. 


5 Event selection 

Events are required to have at least two electrons or two muons within the geometrical accep¬ 
tance regions, and to satisfy the reconstruction, identification, and isolation requirements. The 
pj threshold is set to 20 GeV for the lepton with highest pj, and to 10 GeV for the lepton with 
next-highest pj. The Z boson candidate is formed from the two highest-px^ opposite-charge, 
same-flavor leptons, and must have an invariant mass larger than 50 GeV. In addition, at least 
two jets are required with pj > 20 GeV, within \t]\ < 2.4, and with angular separation relative 
to each lepton of AR > 0.5. The h boson candidate is formed from the two jets that have the 
highest values of the b tagging discriminant. Additional jets in the event are ignored. 
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5 Event selection 


Signal events are expected to have at least one jet b-tagged with the tight and one with the loose 
CSV working points; after b tagging application, the contribution from pileup jets becomes 
negligible. The invariant masses of the two leptons and of the two jets are required to be 
compatible with those of the Z boson (75 <mu< 105 GeV) and h boson (90 < »Jbb < 140 GeV), 
and the value of the £3^®® in the event has to be compatible with zero. After these selections, 
the signal efficiencies range from 10.0% for nip^ = 250 GeV to 20.7% for mp^ = 600 GeV. 

Dedicated control regions are defined to check both the normalizations and distributions of 
the most important backgrounds by inverting the selections used to enhance signal. Drell-Yan 
backgrounds (from qq —> Z/ 7 * —> production) are considered separately as a function 

of the number of b jets, distinguishing Z+jets (no b jets), Z + b (1 b jet) and Z + bb (2 b jets). 
The corresponding control regions are selected by requiring 80 < niu < 100 GeV, £^^®® < 
40 GeV, vetoing of dijet masses close to the Higgs boson mass of 90 < nijj < 140 GeV, and 
applying different b-tagging selections. In the Z+jets control region, no cutoffs are applied on 
b tagging discriminators. The Z + b control region contains events with one b jet fulfilling the 
tight GSV working point, and no other jets passing the loose threshold. Events that enter the 
Z + bb control region should have at least two b-tagged jets, with one passing the tight and the 
other the loose working point. The tt control region is defined by inverting the rtin and £^®® 
selections, dropping the requirement of the dijet mass, and requiring at least one tight and one 
loose b-tagged jet. 

The scale factors that are used to correct the normalization of Drell-Yan and tt backgrounds, 
reported in Table are obtained from a simultaneous likelihood fit to data and simulation 
in the four control regions and are applied in the following steps of the analysis. Multijet 
contamination in control and signal regions is evaluated with data by inverting lepton isolation 
criteria, and its contribution is found to be negligible. The yield of the Vh, ttV, and single top 
production through the s channel are calculated at NLO ||38l EH/ while the other single top 
channels and dibosons are normalized to the measured cross sections ||4m - j42]| . A mismodeling 
in simulation, observed in the control regions relative to data, is corrected by reweightmg with 
a linear function the event centrality distribution (defined as the ratio of the sums in scalar px 
and the energy of the two leptons and two jets in the rest frame of the four objects) in all Monte 
Garlo events, improving the overall agreement between simulation and data. 

Table 1: Scale factors for the four main backgrounds obtained from a fit to the control regions. 
Reported uncertainties are statistical in nature. 


Background 

Z -|- jets 

Z±b 

Z ± bb tt 

Scale factor 

1.069 ±0.002 

0.945 ± 0.012 

1.008 ±0.020 0.984 ±0.010 


An important feature of the signal is that the two b jets originate from the decay of the h bo¬ 
son, whose mass is known with better precision than that provided by the bb invariant mass 
resolution f4^ . The measured jet pj, rj, and (p values are therefore varied according to their 
resolution, in a kinematic fit based on Lagrange multipliers, to constrain the dijet invariant 
mass to wJh = 125 GeV. The fit is used in subsequent steps of the analysis as a discriminant 
in place of mbb- The kinematic fit improves the relative four-body invariant mass resolution 
from 6.3% to 1.2% and 4.0% to 1.9%, respectively, for the smallest and largest values of nip^, 
centering the peaks around their nominal values, as shown in Fig. The effect of the kine¬ 
matic fit is larger at low where the constraint on the Higgs invariant mass has the largest 
contribution. Although both the background and signal m^^b distributions are modified by the 
kinematic fit, the signal significance in a mass window close to the investigated A boson mass 
increases by a factor of two at the lowest mass and by 34% at the highest mass. The resulting 
jet three-momenta are used to redefine all the kinematic variables in the event. 
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Figure 1: Simulated distributions for before (dotted lines) and after the kinematic fits 
(solid lines). Histograms are normalized to unit area. 

Discrimination of signal from backgrounds is achieved through a multivariate discriminant. 
Simulated mass points are divided into three mass regions: low (m^ = 225, 250, and 275 GeV), 
intermediate (wja = 300,325, and 350 GeV), and high mass {niA = 400,500, and 600 GeV). Three 
boosted decision trees (BDT) Il44l are trained separately, one for each region. The inputs of each 
BDT consist of 16 discriminating variables, selected from a list of more than 40 variables: the 
Pj of the Z and h boson candidates, the of tho kinematic fit, the significance of Il36l . 
the dilepton invariant mass, the AR separation and the "twist" angle (defined as tan^^ Acp/Arj) 
between the two b jets ||45l , their GSV discriminator values, the flight directions of the Z boson 
and of the beam in the rest frame of the A boson (cos 6*), the decay angle of the Z boson relative 
to its flight direction in the rest frame of the Z boson (cos 0i), which is sensitive to the transverse 
polarization of the Z boson along its flight direction, the angle of the pull vector Il46l l47l of the 
highest-px jet, which exploits the color connection between the two b quarks originating from 
the h boson, the scalar sum of and the pj of jets and leptons in the event (Sj), the number 
of jets with Pj > 20 GeV, the event centrality, and aplanarity The distribution in BDT 
outputs for data, for simulated signal (S), and for the expected SM background (B) events are 
shown in Fig.|^ weighting each entry in the m££\,j distribution by the expected S/(S+B) ratio 
of the BDT bin in the signal-sensitive region with BDT > 0.6. 


6 Systematic uncertainties 

The sensitivity of the analysis is currently limited by the available data, and not by systematic 
uncertainties. 

The uncertainties in the normalizations of the four main backgrounds (Z+jets, Z + b, Z -|- bb, 
and tt) originate both from the fits in the control regions and from the extrapolation to the signal 
region. The former are reported in Table and the latter is evaluated through a simultaneous 
likelihood fit to data and to simulated yields in several statistically independent regions, ob¬ 
tained by altering the selections used to define the four control regions. A 13% normalization 
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6 Systematic uncertainties 



Low mass BDT 
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Figure 2: BDT outputs and invariant mass distributions in the low, intermediate, and high mass 
regions. The plots are for BDT > 0.6, weighted by S/ (S + B) in each BDT bin. Histograms 
for signal are normalized to the expected exclusion limit at 95% confidence level. Statistical and 
systematic uncertainties in simulated samples are shown as well. Either the ratio (left) or the 
difference (right) between data and SM background is given at the bottom of each panel. 
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uncertainty due to this extrapolation is assigned to Z+jets, 12% to Z + b, 2.1% to Z + bb, and 
6.2% to tt events. Normalization uncertainties for other SM backgrounds correspond to the 
ones on their measured or theoretical cross sections. 

The uncertainties in lepton reconstruction, identification, isolation, and trigger efficiencies are 
evaluated through specific studies of events with masses in the region of the Z peak. Uncertain¬ 
ties in background or signal normalization and distribution due to uncertainties in jet energy 
scale and resolution i3? ll and b tagging scale factors 11371 are estimated by changing the corre¬ 
sponding values by ±1 standard deviation (cr). Additional systematic uncertainties affecting 
the normalization of backgrounds and signal from the choice of PDF Il48l |4^ , contributions 
from pileup, fluctuations because of the presence of unclustered energy in the event, and 
integrated luminosity |T7 i| are also considered in the analysis and are reported in Table [^(Nor¬ 
malization). 

Results are extracted from an analysis based on a binned likelihood fit to the two-dimensional 
(2D) distribution of versus BDT. Dependence on jet energy scale and resolution, b tag¬ 
ging, as well as factorization and renormalization scales are propagated to the 2D-templates, 
taking into account the correlations between the two variables. The impact of reweightmg the 
Monte Carlo distributions is considered as an additional source of background in modeling the 
uncertainty. Finally, the uncertainty from the limited number of simulated events is treated as 
in Ref. liSOl . The sources of systematic uncertainty affecting the forms of the distributions are 
summarized in the first column of Table [^ (Shape). The second and third columns indicate the 
respective ranges in the relative impact in percent, obtained by the changes implemented in the 
background and signal contributions. 

The systematic uncertainty with the largest impact on the expected limit is from the reweight¬ 
ing of the background (accounting for a <6% difference on the expected limit, depending on 
niA), which is followed by the limitations in the number of simulated events (^4%), and the 
factorization and renormalization scales (<2%). The effect of the other sources is small (^1%). 


7 Results and their interpretation 

Results are obtained from the combined signal and background fit to the binned two-dimensional 
distribution of the four-body invariant mass arid the BDT output in the signal-sensitive 
(BDT >0.6) region. With no evidence of significant deviation from background expectations, 
the asymptotic modified frequentist method is used to determine the limit at the 95% confi¬ 
dence level (CL) on the contribution from signal, treating systematic uncertainties as nuisance 
parameters that are integrated over in the fit ll5T] - l54l . 

The observed limit, as well as the expected limit and its relative ±1 and ±2cr bands of uncer¬ 
tainty, are reported as a function of the A boson mass in Fig.j^for cta f3(A ^ Zh —> ££bb), i.e. 
the product of the cross section and the A —> Zh, h —> bb, and Z ^ branching fractions, with 
£ = e or ji. The limits are obtained by considering the A boson produced via the gluon-gluon 
fusion process in the narrow-width approximation. Interpolated mass points are obtained as 
in Ref. Il55l , and numerical values are reported in Table [^ A signal upper limit at 95% CL is 
set on aA B{A —> Zh —> ££hh), excluding from 10 to 30 fb for niA near the kinematic threshold, 
~8 fb for niA ~ 2mt, and up to ~3 fb at the high end (600 GeV) of the considered mass range. 
Comparable limits have been recently obtained for the same channel by the ATLAS Collabora¬ 
tion llT4ll . The most significant excess at m^^bb = 560 GeV has a local significance of 2.6c7', which 
becomes a global l.ltr significance taking into account the look-elseivhere effect 11561 . 
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7 Results and their interpretation 


Table 2: Summary of systematic uncertainties. Normalization: sources of systematic uncer¬ 
tainty and their effect in percent on the normalization of signal and background distributions. 
Shape: sources of systematic uncertainty and the range of their effect in percent on relative 
changes made in the form of the background and signal distributions. 


Sources 

Backgrounds 
Drell-Yan, tf Others 

Signal 

Normalization 

Gontrol region fitting 

<2.4% 

— 

— 

Extrapolation 

2-13% 

— 

— 

Lepton and trigger 


2 5% 

2 5% 

efficiency 




Jet energy scale 

— 

5.7% 

3.8-0.2% 

Jet energy resolution 

— 

3.2% 

0.8-0.5% 

b tagging 

— 

4.9% 

3.6-3.2% 

Unclustered energy 

— 

1.9% 

1.4-1.0% 

Pileup 

— 

0.9% 

1.2% 

PDF 

— 

4.3% 

4.0-7.9% 

Gross section 

— 

9.2-15% 

— 

Integrated luminosity 

— 

2.6% 

2.6% 

Shape 

Jet energy scale 

<4% 


<8% 

Jet energy resolution 

<2% 


<4% 

b tagging 

<4% 


<8% 

Factorization and 




renormalization scales 

<6% 


6-10% 

Monte Garlo reweighting 

<15% 


— 

Monte Garlo statistics 

1M% 


— 


For wja > 2mt, the width of the A boson depends strongly on the model parameters. Different 
limits are provided by taking into account the natural width of the A boson (Fa) in the recon¬ 
structed leaving the BDT unchanged. Figure]^ shows the exclusion limit above 2m t for 

an average width of 30 GeV, and the dependence of the observed limit on Fa for wja = 500 GeV. 
The local and global significance at = 560 GeV become, respectively, 2.9(7 and 1.5cr, as¬ 

suming Fa = 30 GeV. 

As an independent cross-check of the 2D fit, the signal is extracted by applying two comple¬ 
mentary strategies, based on one-dimensional fits. The first one consists of fitting the 
distribution, after selecting events in a signal-enriched region by applying a BDT > 0.8 se¬ 
lection. The second relies on fits to the BDT distributions, after selecting events within the 
resolution of the signal m^^bb peak. The two methods give upper limits compatible with those 
from the 2D fit, but 10 to 20% less stringent. 

The results are interpreted in terms of Type-I and Type-II 2HDM formulations i4J. The S(A —> 
Zh) and 6(h —> bb) branching fractions and the signal cross sections are computed at next-to- 
next-to-leading-order (NNLO) with SuShi 1.2.0 157 1 and 2hdmc 1.6.4 Il58l . respectively, using 
the MSTW2008LO, NLO, and NNLO sets of PDF. The S(Z —)■ ££) branching fraction, with £ = e 
or is taken from the measured value 1591 . Both gluon-gluon fusion and associated production 
with b quarks have been considered. The latter is rescaled to the fusion process, taking account 
of the difference in acceptance for signal, as well as the efficiency for selecting dijet pairs in the 
presence of combinatorial contributions from additional b quarks in the event. The parameters 
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used for the models are: mh = 125GeV, ntu = ?nH± = wa, m ^2 ~ +tan^/l], 

^ 6,7 = 0 IMl/ while niA = 225-600 GeV, 0.1 < tan/I < 100, and —1 < cos(j6 — a) < 1, using 
the convention 0 < fi — a < 7i, where tan fi and a are, respectively, the ratio of the vacuum 
expectation values, and the mixing angle of the two Higgs doublets IH. 


The observed limit, together with the expected limit and its relative ±1 and ±2a uncertainty 
bands, are shown in Fig. interpreted in Type-I and Type-II 2HDM for niA = 300 GeV. A 
sizeable fraction of the 2HDM phase space is excluded at a 95% GL with respect to the previous 
GMS searches 1131 . 


Table 3: Observed and expected 95% GL upper limits on ox 'S(A —> Zh ^ f'dbb) as a function 
of niA in the narrow-width approximation, including statistical and systematic uncertainties. 


rriA [GeV] 

225 

250 

275 

300 

325 

350 

400 

500 

600 

Observed [fb] 

17.9 

16.8 

14.8 

19.5 

10.1 

8.84 

3.29 

3.35 

2.61 

Expected [fb] 

17.9 

18.1 

16.4 

13.6 

10.0 

7.84 

5.27 

2.79 

1.93 

A-> 

Zh-> 

eebb 




L = 

19.7 fb- 

(8 TeV) 


95% CL limits 

-Observed 

. Expected 

Expected ± 1a 
Expected ± 2a 
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450 


500 
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[GeV] 


Figure 3: Observed and expected 95% GL upper limit on ux S(A —> Zh ^ ££bb) as a function 
of WJA in the narrow-width approximation, including all statistical and systematic uncertainties. 
The green and yellow bands are the ±1 and ±2(7 uncertainty bands on the expected limit. 


8 Summary 

A search is presented for new physics in the extended Higgs sector, in signatures expected from 
decays of a pseudoscalar boson A into a Z boson and an SM-like h boson, with the Z boson 
decaying into (E being either e or /r) and the h boson into bb. Different techniques are 
employed to increase the sensitivity to signal, exploiting the presence of the three resonances A, 
Z, and h to discriminate against standard model backgrounds. Upper limits at a 95% GL are set 
on the product of a narrow pseudoscalar boson cross section and branching fraction o'aB{A ^ 
Zh —)• EEbb), which exclude 30 to 3 fb at the low and high ends of the 250-600 GeV mass range. 
Results are also presented as a function of the width of the A boson. Interpretations are given 
in the context of Type-I and Type-II 2HDM formulations, thereby reducing the parameter space 
for extensions of the standard model. 
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8 Summary 
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Figure 4: Observed and expected 95% CL upper limit on o'aB{A —> Zh —> £tbb) for Fa = 
30 GeV as a function of niA (left), and for niA = 500 GeV as a function of the width of the A 
boson (right). 
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Figure 5: Observed and expected (together with ±1, 2a uncertainty bands) exclusion limit for 
Type-I (left) and Type-II (right) models, as a function of tan/l and cos(/3 — a). Gontours are 
derived from the projection on the 2HDM parameter space for the niA = 300 GeV signal hy¬ 
pothesis; the observed limit is close to la above the expected limit, as shown in Fig.|^ 
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